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Abstract 

We calculate expectation values of QCD operators consisting of the products of 
the four operators of the light quarks q~T x qq~T Y q, with T X ' Y corresponding to the 
scalar, pseudoscalar, vector, pseudovector (axial) and tensor Lorentz structures, in 
the nucleon. All combinations of the light flavors are considered. For the evaluation 
we use elements of the Perturbative Chiral Quark Model (PCQM), approximating 
the contribution of the valence quarks by the contribution of the PCQM constituent 
quarks. The contribution of the sea quarks is treated by averaging over the QCD 
pions with the distribution of the pion field being determined by the PCQM. For 
quarks with the same flavor the expectation values are dominated by the contri- 
bution of the sea-quarks. In the scalar case the contribution of the sea quarks is 
dominated by the "disconnected terms" were one of the pairs of the qq operators 
acts on the vacuum while the other one acts on the quarks of the pion. The role 
of the interference terms with one of the qq pairs acting on the sea quarks and 
another one acting on the valence quarks increases for quarks with different flavors. 
The result for the scalar condensate is compared to the one obtained earlier in the 
framework of the Nambu -Jona-Lasinio model. 



1 Introduction 

It is known that expectation values of the two-quark QCD operators qq give the total 
number of quarks and antiquarks in a hadron under certain reasonable assumptions jTJ 
12]. The motivation for studies of the expectation values of the four-quark operators in 
hadrons is that they carry information about the correlation of qq pairs. These expectation 
values in nucleons determine the coefficients of the next-to-leading order operators of deep 
inelastic scattering [3]. One more reason is the manifestation of such operators in QCD 
sum rules in nuclear matter jl] |3] . The lack of data on the expectation values of these 
operators became one of the main obstacles for further development of this approach. 
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Unfortunately, the calculations of these expectation values require some model as- 
sumptions on the quark structure of the nucleon. As it stands now, the only calculation 
of the four-quark condensates in the nucleon is the one carried out by Celenza et al. for 
the scalar case in the framework of the Nambu and Jona-Lasinio model j7| under certain 
additional assumptions. Also, some constraints on the values were obtained by Johnson 
and Kisslinger jS] by analyzing QCD sum rules for nucleons and isobars. 

In this paper we calculate the expectation values of the four-quark condensates in 
nucleons by using elements of the perturbative chiral quark model (PCQM). The chiral 
quark model, originally suggested in P, was fully set up in (TUJHUII2]- In the PCQM 
the nucleon is treated as a system of relativistic valence quarks moving in an effective 
static field. In addition, the valence quarks are supplemented by a perturbative cloud 
of pseudoscalar mesons as dictated by chiral symmetry requirements. In the current 
paper we restrict to the simplest SU(2) version of the PCQM, which includes only pions. 
To facilitate the evaluation of the four-quark operators we resort to previously derived 
results [IT] , such as nucleon wave function renormalization and self-energy contributions, 
which are used as an input in the present derivation. Although we do not apply the full 
perturbative machinery of the PCQM, as laid out in ^0], the present evaluation serves as 
a first indication for the values of the four-quark condensates. 

We calculate the expectation values of the four-quark operators as the matrix elements 
of the PCQM nucleon. In general the QCD operators q act on the valence and the 
sea quarks. We approximate the averaging over the valence quarks by averaging over 
the PCQM constituent quarks and also approximate the averaging over the sea quarks 
by averaging over the pions. In the original formulation of the PCQM the pions were 
treated as separate degrees of freedom, without taking into account their quark structure. 
Actually, in ^TJ the sigma-term was calculated in the framework of PCQM, with the 
pions determining the sea-quark contribution. We shall employ an additional extension 
of the PCQM, considering the pions as physical particles with their quark contents being 
determined by QCD. In this approach we therefore obtain the excess of the four-quark 
operator expectation value over the one of the vacuum in the nucleon volume. 

In the following we calculate the expectation values 

v XY,hf 2 = ( N \ T XY ^f*\N) (1) 
of the four-quark operators of light quarks 

t xy ^ = (: g -A<rV ia ' • W 0(<W<W - <WM- (2) 

Here T x , T Y are the matrices acting on the Lorentz indices, q* denote QCD quark opera- 
tors with / standing for the flavor. The dots denote the normal ordering of the operators, 
a, a', b, b' represent the color indices. We calculate the condensates for the basic 4x4 
matrices 

V s = /, r p * = 75 , ^ = 7m , T^ = 7m75, r p =^(7 M 7,-7,7,) (3) 
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with 7^ being the Dirac matrices, that is we consider the scalar, pseudoscalar, vector, 
pseudovector (axial) and tensor cases. We also obtain the results for the mixed condensate 
jjsv,du which i s important in applications. 

These are three types of contributions to JJ XY ^ lf2 in our approach. All four operators 
q can act on the constituent quarks, providing the term C xy '^ 2 . Also, four operators can 
act on the pions providing the term p^/ 1 / 2 . There is also a possibility that two of the 
operators act on the constituent quarks while the other two act on the pions. Denoting 
the last term as J XY ^f 2 , we present the expectation values as 

JjXY,hf 2 = QXY,hh + pXY,hf 2 + jXYJih. (4) 

In Appendix A we show how these contributions manifest themselves in the PCQM for- 
malism. 

To simplify the notations we introduce 

with a similar convention for the other functions (C, P, J, T) involved. 

While the four-quark condensates are Lorentz scalars in the scalar and pseudoscalar 
channels, they have a more complicated structure in the case of the vector and axial 
channels: 

U^ = a v ^ + b v ^^. (6) 

Here is the momentum of the nucleon, while m denotes the nucleon mass. Also, in the 
tensor channel we have 

= a s ^,a/3 + b t^ a/ 3 (7) 

with 

Sfiu,a/3 = ~ 9^/39ua, (8) 

V,"/3 = + PvPpg^a - P^Ppdva ~ PuPaQ^)- 

We shall denote the values of a v ' A ' T and b v ' A ' T corresponding to the contributions C, P 
and J by the lower indices, i.e. a^p j, etc. 

We approximate the averaging of the operator qT x qqT Y q over the valence quarks by 
the expectation values of the products of the constituent quark operators Q averaged 
over the renormalized constituent quark PCQM states. Due to the normal ordering of the 
operators they should act on two different quarks. The expectation value is proportional 
to the probability to find the two constituent quarks at the same space point 



/ mr)\ 4 d 3 r 



AttR 3 



with R standing for the size of the system of the three constituent quarks (quark core 
radius), while ip is the constituent quark wave function. The factor comes from the 
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four angular wave functions 7^172 integrated over the solid angle. Since the PCQM deals 
with a quark wave function ip(r) provided in explicit form, the contributions (7 >-' 1 -' 2 are 
also evaluated explicitly. The effect of the wave function renormalization induced by the 
interaction of the constituent quarks with the pion provides noticable corrections in the 
case of scalar and vector structures only. 

The averaging of the operators qT x qqT Y q over the sea quarks is treated as the expec- 
tation value of these operators in pions. The distribution of the pion field is assumed to 
be the one determined by the PCQM. 

The pion expectation values where expressed in paper [T3] by using the current al- 
gebra technique through the four-quark expectation values in vacuum. We obtain the 
expressions for P XY >hfe through these expectation values. However, to obtain the spe- 
cific numbers, we use the factorization approximation for the vacuum expectation values, 
suggested first by Shifman et al. [Tl]. In the factorization approximation it is assumed 
that the vacuum states dominate in the sum over the intermediate states. While there 
are indications that this approximation may be violated in some of the channels [3], the 
factorization was advocated recently in jT^j. Under this approximation the expectation 
values P XY -^h are expressed by the qq vacuum expectation values which are known to 
be PS] 

M 2 F 2 

qq = (9) 

m u + m d 

for each of the light flavors with M n and F n being the mass and the decay constant of the 
pion, while m u ^) are the current quark masses. Here we adopt the notations accepted in 
chiral perturbation theory [T2j. In Eq.Q q stands for the u or d quark field and isotopic 
invariance of the vacuum is assumed. The value (0|gg|0) gives the characteristic size of the 
contribution of the pion sea. Thus the contribution of the constituent quarks is expected 
to be smaller than that of the sea quarks since ^| (0|gg|0)i? 3 | ~ |. 

In the "interference term" one of the operators qT x q acts on the valence quarks while 
another one qT Y q acts on the sea quarks. Following our strategy, we approximate the 
corresponding matrix elements by those averaged over the constituent quarks and over 
the PCQM pion field. There are several possibilities to insert this four-quark operator. 
The operator can connect the pion with any of the constituent quarks of the nucleon. 
This contribution (the "contact interference") is proportional to (7r|gr y g|7r) which does 
not vanish for the scalar case only. Thus, only the contributions J sv = J vs and J have 
nonzero values which are proportional to the expectation value |T7j 

Wfthr > = = _M. (10) 

However, there is an additional small factor besides the characteristic parameter (0|gg|0). 
This factor reflects the small probability for the pion and the constituent quark to overlap 
at the same space point for a nucleon Fock state described by the valence quarks and a 
pion. The interference process can determine a QQn vertex as well, since the pseudovector 
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and pseudoscalar currents connect the pion and vacuum states. Thus, one can consider 
the self-energy diagram with one of the PCQM vertices being replaced by the four-quark 
operator. That would be the first order diagram in the PCQM 7rQ interaction. This 
"vertex interference" leads to a numerically larger contribution, except for the case of 
the scalar-vector condensate. The neutral pions provide the contribution to the terms 
J A, ff and J Ps 'ff of the quarks with the same flavor. The charged pions contribute to the 
condensates uY ddY u and uY Ps ddY Ps u. Thus, they provide the contributions to all the 
structures J x > ud with coefficients defined by the Fierz transform. 

Following the general strategy of the PCQM we include only the lowest order con- 
tributions in the irQ interactions. We also assume that only the ground states of the 
constituent quarks are included as intermediate states in the self-energy diagrams. This 
means that the nucleon and delta-isobars only are included as intermediate states of the 
nucleon self-energy. This is a standard assumption of PCQM calculations [PUl lllj. 

It was shown in ^Hj that in the scalar case in the four-quark pion expectation value 
a "disconnected term", in which one of qq pairs acts on vacuum, can be singled-out in 
a natural way. This is strongly pronounced in the case of the color-singlet four-quark 
operator q a q a q b q b for which 

(ir\qqqq\ir) = 2(0|gg|0)(7r|gg|7r) + (n\(qqqq) int \7r} . (11) 

The "internal" contribution presented by the second term on the right hand side (rhs) 
of Eq. lfTTj) appeared to be about an order of magnitude smaller then the "disconnected" 
one, presented by the first term. For the color asymmetric operator determined by Eq. (J2J) 
the "disconnected" terms still provide about | of the total pion expectation value. This 
leads to the natural presentation. 

p SJlh = pSd lh + p?J lh (12) 

with 

PS^=l-2{0\gg\0){N\(qq) Ma \N) (13) 

for the different flavors, where \N) represents the fully dressed nucleon state. The factor 
| (it is | for identical flavors) on the right hand side of Eq. (|13j) is the weight of the 
colorless combination q a q a in the color asymmetric expectation value defined by Eq.(J2J. 
In the case of the scalar condensate the sea quarks provide the main contribution. A large 
part of it is determined by the disconnected term, related by Eq. (jl3|) to the sea-quark 
contribution to the nN sigma term. The internal contributions, coming mostly from the 
sea quarks are several times smaller. The expectation values of the operators with the 
same flavor (uYu) 2 in the other channels are also determined mostly by the sea-quarks. 
For the mixed-flavor condensate uYudYd the sea quark terms and the interference terms 
provide contributions of the same magnitude in most of the channels. The valence quarks 
provide a smaller correction. In the case of the scalar-vector condensate P = 0, and 
for the neutron the valence quarks provide the main contribution, while in the proton the 
interference effects contribute of the same order. 
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In the calculations carried out below, we use the values F w = 92.4 MeV for the pion 
decay constant and the value m u + = 11 MeV for the sum of the light quark masses. 
Latter value, given in [THj, leads to the conventional value (0|gg|0) = (—245 MeV) 3 at the 
normalization scale of 1 GeV. This set of values was also used in the paper [U] . Note that 
in the papers fDJ another value for the sum of the quark masses has been used, e.g. 
m u + rrid = 14MeV. This value was also given in [TSj as one of the possible ones. Both 
values for m u + are consistent with nowadays experimental data [T9"j . 

We present the results for the condensates (uT x u) 2 both for proton and neutron. The 
values of the (dT x d) 2 condensates are determined by the isotopic invariance relations 

(p\(dT x d) 2 \p) = (n\(uT x u) 2 \n), 
(n\(dT x d) 2 \n) = (p\(uT x u) 2 \p) 
while for the mixed flavor condensates we have 

(p\uT x u ■ dY x d\p) = (n\uT x u ■ dY x d\n) 

except for the pseudoscalar case Y x = T Y = 75. In the pseudoscalar case an explicit de- 
pendence of the interference terms on the current quark masses m Uj d causes contributions 
which break the isotopic invariance. These terms are numerically small. 

The results enable to obtain also values for the condensates uT x d ■ dT x u. This can be 
done by using the Fierz transform. 

We compare the value of the contribution (N\ (uu + dd) 2 \N) with the value obtained in 
jH] in the framework of the NJL model under certain additional assumptions. The values 
appear to differ by about 70%. 

For the sake of simplicity we shall use the wording "scalar", "pseudoscalar", etc. con- 
densates for the expectation values of the operators with the repeated Lorentz structures 
qTqqVq. Thus the "scalar" expectation values are rather "scalar-scalar" ones, etc. 

2 Contribution of the valence quarks 

In this Section we calculate the contribution to the four-quark expectation values arising 
from averaging over the system of three valence quarks. Using the results of Appendix A 
we present this contribution by 

c XYj lh = f ^xq h {x)Y x q h {x)q h {x)Y Y q h {x)\(j) ) (14) 

with 0o denoting the nucleon as a bound state of three valence quarks. Our main assump- 
tion here is that the matrix element in the rhs of Eq. (jl4j) is approximated by the matrix 
element of the renormalized constituent quark operators Q r , i.e. 

c xY Jlf2 = J d 3 xQ r > fl {x)r x Q r > fl {x)Q r > h {x)r Y Q r > f2 {x)\<f> ) (15) 
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The renormalization effects are expected to manifest themselves through small cor- 
rections only. As we shall see below, these corrections are of the order of several percent 
only, except for the scalar and vector structures. Thus, we start with the unrenormalized 
constituent quark operators Q in which the pion cloud is not included. The corresponding 
contribution 

C xyjl/2 = (M J d 3 xQ h (x)T x Q^(x)Q h (x)T Y Q^(x)\ ( l ) o) (16) 

is illustrated by Fig.l. The constituent quark operators provide a nonzero value while 
acting on different quarks of the 0o system only. This is due to their normal ordering. 
Thus we find immediately 

Cn Y ' W = (17) 

for the neutron. Using Eq. (jl6j) we obtain expressions for the contributions of the con- 
stituent quarks through the wave functions ipi(x). Assuming that the constituent quarks 
U and D are described by the same wave functions ijj u {x) = ipd(x) = we present the 
general expressions for the proton as 

Cf uu = J d 3 xF(x) (18) 

with 

F(x) = Pi){x)T x i){x) ■ i){x)T x i){x). (19) 

Here P stands for the projection on the symmetric spin state of the two-quark system, 
while the total antisymmetrization is provided by the color variables. For the condensate 
uT x udT x d we find both for proton and neutron 

C x ' ud = 2 J d 3 xJ r (x) (20) 

since there are two ud pairs. 

The invariant coefficients of the rhs of Eqs.(JB} and (|IJ) can be obtained in a specific 
reference frame. Assuming i and j to be the three-dimensional indices, corresponding to 
the four-dimensional indices fi and z/, we find in the rest frame of the nucleon 

a c = -\c* X ' hh 5 l3) b c = C xx ' hh -a c (21) 

for the coefficients of Eq.©, i.e. for the vector and pseudovector cases. Denoting the 
three-dimensional indices corresponding to the four-dimensional indices a and j3a.sk and 
/, respectively, we obtain in the same frame 

a c - a^ijkl d ik°jl, O c - -a c - -OjlL'Qj^ . 
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Above equations are true for any constituent quark model. In the specific case of the 
PCQM the wave functions of both U and D constituent quarks are [TT] 

= WeC "*' ( ipkx ) (23) 

with the normalization constant 

N= [tt 3 / 2 ^ 3 (1 + |/5 2 )]~ 1/2 (24) 

and x being the two-component spinor. 
The model parameters 

/3 = 0.39, R = (0.6 ± 0.05)fm (25) 

are fitted to reproduce the value of the axial coupling constant and of the proton charge 
radius. We will present the numerical values for the mean value of R=0.6 fm. 

o 

A straightforward calculation provides for the expectation values Cp " in the proton 



3 „ . 15 

with 



Cf"" = (1 - + (26) 



while the value is zero for neutron. For the uudd condensate we get 

ftS,«d = 2(1 _ ^2 + 1^4^2_ (2g) 

2 16 

The details of the calculations for the other structures are presented in the Appendix B. 
For the structures (uT x u) 2 in the proton we find for the pseudoscalar case 



'v 

while in the vector channel 



c Ps,uu = _( 3 2 U 2^ ( 29 ) 



= -^ 2 > b v c , p = (1 + y/? 2 + ^/3 4 )AT 2 , (30) 
and for the pseudovector case 

al P = -\{l-\p+ l ^)N\ k, P = ~°4, P - (31) 
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Note that be = — Oq since the matrix element of the time component of the pseudovector 
operator turns to zero. This is true for the solution of the Dirac equation in any effective 
field. For the tensor case we get 

al = \{l + \f+ l ^)M\ bl = -\{l + \f+ l ^)M\ (32) 

Following the previous analysis these values turn to zero for the neutron. 

Turning to the case of different flavors, we find the expectation values of the operators 
uY x udY x d in a nucleon to be twice as large as the values of {uY x u) 2 in the proton 

C x > ud = 2 Cp' uu . (33) 

We also present an example of the condensate for the mixed scalar-vector structure 
rpsv,du _ ddu^^u, which is needed in applications. In the rest frame of the nucleon only 
the time component of the vector T M survives, providing 

C sv ' du = 2(1 - ^/? 4 )AT 2 . (34) 
16 

It is convenient to express the values in "units" of the value Eq, defined as 

£ 3 = -(0\qq\0), e ~ 245 MeV. (35) 

To get a feeling for the relative size of the contributions, we present the numerical 
value 

N" 2 = 1.50 • l(T 3 GeV 3 ~ O.lOe 3 (36) 

which is the result of the straightforward computation of the rhs of Eq. (j27|) . 

The 7rQ interactions provide the nonzero values of the condensates uY x uuY Y u in the 
neutron. This happens since the four-quark operators connect the only valence U quark 
with the intermediate U quark of the n~U self-energy loop of the valence D quark. The 
value is 

SlX.UU O /■iX.UU 

n ~ ~ d¥ v ^ ' 

with S~ standing for the contribution of 7r~ to the self-energy of the valence quark with 
the energy E. The direct calculation gives = —0.082. 

Now we take into account the change 5ip of the shape of the single quark wave function 
ip(x), caused by renormalization [TO] 



ip r (x) = tp(x) + 5tfj(x), 5i/j(x) = A(h(x) + 7o)^(jc) (31 

with 

_ 5m (3R | + f/3 2 x 2 

2 l + f/3 2 ' m ~ l + l(3 2 R 2 ' 



Here 5m < is the shift of the effective mass of the constituent quark caused by the 
pion cloud. The numerical values are A = —2.5 • 1CT 2 , h(x) = 1.06 — j^-. 

The term containing the function h(x) provides corrections which do not exceed 3%. 
This happens due to the strong cancellations of the two terms, composing h(x). We shall 
neglect these corrections. The term containing the Dirac matrix 7 mixes the Lorentz 
structure of the condensates. It provides nonvanishing contributions to the scalar, vector 
and scalar-vector expectation values. 

Thus, we obtain for renormalized values defined by Eq. (|15jl 

C^ uu =C% UU + 4A Cl v ' uu - , C s ' ud =C s ' ud + 2A(C sv ' du + C sv ' ud )- , (39) 

P P P ^ ry 

c sv,du w gsv,du + 2K(c SM l+ be -) 

7 

with 7 = Y+i/P"- Thus, the scalar-scalar condensates are reduced by 16% due to the 

renormalization effects. The scalar-vector condensate is reduced by 11%. Also, in the 
vector case we have 

«c=«c, 6g=6S(l + 4A 7 ) (40) 

o 

reducing the value of b c by about 6%. For the pseudoscalar, axial and tensor structures 
the corrections are negligibly small and we put 

in these cases. 



3 Contribution of the sea quarks 

Now we calculate the contribution of the sea quarks. In the PCQM the excess of the sea 
quarks in nucleons over the QCD vacuum sea inside the nucleon volume is contained in 
the mesons, coupling to the constituent quarks. In the SU(2) version of the model, which 
we assume in this paper, only the pions are included. In the framework of the PCQM 
this contribution is contained in the next-to-leading order of the model. In other words, 
it is sufficient to include pion exchange in the one-loop approximation. 

The distribution of the pion field 7r a (x) is determined by the PCQM quark-pion in- 
teraction 

Hj(x) = ^(x)ij 5 Six)T ^ a{x) ^(x) (42) 

where Hf(x) represents the SU(2) doublet of light quarks, while S(x) is the effective scalar 
field. 
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In the one-loop approximation of the PCQM the pions are contained in the constituent 
quark self-energy diagrams and in the diagram, describing the pion exchange between 
the constituent quarks (Fig. 2). The contribution of the sea quarks (see Eq.(|TJ)) can be 
presented as 

^'-d.t™io(| + |) (43) 

with S a and A" standing for the self-energy and exchange contributions of the pions ir a 
(a = +, — ,0). A similar presentation was actually used in [TT] for the calculation of the 
sigma term. In that case T was the scalar quark operator qq. 



The rhs of Eq.(j43j) can be simplified by noticing that in the PCQM the relation for 
the total energy shifts caused by the self-energy and exchange diagrams 20J 

EA a = ^E^ (44) 

a ^ a 

holds also for each pion a separately, when limiting single quark lines to the ground state. 
Using Eq. (j44|) we present the total pion contribution to the nucleon mass as 

10, 



This leads to 



for the proton with 



p xYj lh = 3/3° + 4/3+ + 2/T ggj 



9 dMl (45) 



(3 a = (7i a \T XYJlh \7i a ). (46) 



The coefficients multiplied by the expectation values (5 a are determined by the numbers 
of the quarks which can emit the pion 7r a and by the strength of the ttQQ vertex. For 
example, 4/3 + means that there are two quarks (these are ?7-quarks) coupling to a n + , 
and each of the n + DU vertices contributes the factor V% etc. 

For the neutron we get 

P xr, hh = 3/3° + 2/3+ + 4/3- gg, 

9 dMl' 1 ] 

The value was evaluated earlier in the calculation of the sigma term ^T], providing 

(W 

« l^eV- (48) 



The pion expectation values (3 a can be expressed by the vacuum expectation values of 
the four-quark operators. This was done in JH] by using the reduction formula obtained 
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by Lehmann, Symanzik and Zimmerman [21] . Due to the partial conservation of the axial 
current (PCAC) the pion state vector can be expressed by the vacuum as (see e.g. |2*2*j ) 

M = 7^^ (I)|0> - (49) 

Here A® 5 (x) is the axial current of the light quarks 

A; i5 (x)=Y.d c (xh, l5 u c (x) (50) 

c 

with c being the color index. (We shall assume the summation over the colors in all the 
equations presented below.) This enables to present the expectation values defined in 
Eq. (|46|) by the vacuum matrix elements [13] 

^ = JL(0|B«|0) (51) 

71 

with 

B a = ^yj d 3 xdy dz 5(x - y )S(z - x Q )[Q%(zo), [Q^y ),T XY ^(x)}}. (52) 

Here V is the normalization volume and in the double commutator occur the axial charges 
Q5, corresponding to the axial current A^ 5 (x). For example, in the scalar channel it was 
found 

B ± = -{u a u b u a 'u b ' + u a u b d a 'd b ' + d a lb u b u a ' lb d b '){5 ah 5 alv - <W<U) (53) 

for the operator T s,uu , averaged over the charged pions ir ± . The total contribution of the 
pion cloud to the uuuu condensate in the SS channel - see Eqs. (J45|) . (JjTj) - takes the form 

P s ' uu = 5L(0|2w a wW + u a u b d a 'd b ' + u a ^u b u b '^u a ' + u a ^ 5 d b d b '^u a ')\0) (54) 

3F W 



dMl 



In ^3] the pion expectation values were expressed by the ones of the vacuum for all 
channels. Thus similar equations can be presented for all structures. Since only the 
vacuum expectation values are involved, we find for the coefficients in the rhs of Eqs.([B|). 





b^' T = 0. (55) 

To avoid complicated formulas we shall present the final results in the framework of the 
factorization hypothesis for the vacuum expectation values. 

It is convenient to present 

2 1 

<W<W — 5ab>5ba' = ^aa'hb' — X! Kia'^bb' (56) 

1 p 
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with A p standing for the SU(3) Gell-Mann matrices normalized by the relation Sp A p A r = 
25 pT . 

In the factorization approximation we find for quarks of the same flavor 

(0\qT r qqT s q\0) = l[S P r r S P r s - is P (r r r s )] ((0|gg|0>) 2 (57) 
Id 6 

for any 4x4 matrices r s r acting on Lorentz indices. If the quarks have different flavors 
we come to 

(o|%r r ^r s g,|o) = ^s P r r s P r.(o|g i9i |o)(o|^|o) (58) 

and 

(o|g i r r ^r aft |o) = -^s P r r r a (o|ft ?i |o)(o|^|o). (59) 

For the matrices T p s = T rjS X p this approximation provides 

(oIE^^t^io) = -is P (r r r a )<o|g <ft |o)(o|^|o> (eo) 
p J 

which is true for i — j and i ^ j, while 

(0|5>ftftfcfZ<&|0> = (61) 
p 

for i 7^ j. 

In the factorization approximation the contribution of the sea quarks contains the 
f cLctor 

(M£ = _^<0 H 0) (62) 
- see Eq.tjUJ), and we can present 

pX,hf 2 _ M l 3 qX,hf 2 



u + m d dMl p ' n 



Here we denoted Pp* = P*^ 1 ? 2 and the rhs of Eq. (|4fi|) turns to zero for 1^7. The 
lower index denotes proton or neutron. 

Using the results of [13], we find for the same flavors 



qS,uu qS,uu qPs,uu qPs,uu ' (P.A\ 

°p — °n — ~g~~5 °p — °n — g5 1011 



2 2 

qV,uu _ qA,uu qV,uu _ cA,uu _ _ 

D p ~ D p ~ g"/^> n ~~ n ~~ gy^f) 

nT,uu oT,uu 

13 



with the tensor s^ a/ 3 defined by Eq.fjHJ). 

For the quarks of different flavors we obtain for the proton and neutron 
14 2 2 2 

qS,ud _ nPs.ud nV,ud _ qA,ud nT,ud fKK\ 

*-> — ~~ "g"> ° — 9' _ ' ~~ ~~ _ g S V,a/3- ^OOJ 

We now show that in the scalar channel the disconnected terms are separated in a 
natural way. We have the result 



qS,uu _ 5 1 Q S,ud _ ' 2 , . 

D p,n — 3 g' _ 3 9' 



As it was shown in [T3] the expectation values of the scalar four-quark operators are 
dominated by the disconnected terms with one of the qq pairs coming from the vacuum. 
This corresponds to the approximation 

k 4 

qS,uu quu qS,ud qud / >y\ 

°p,n ~ °dis ~ g > D p,n — °dis ~ g • K° 1 J 

On the other hand, the factor m M + in the rhs of Eq. fjo^j) is just the expectation value 
of the operator qq in the pion - see Eq. (110)1 . Then we have 

M 2 f)V f)Y 1 

f ' = <*lwk> = ~(N\uu + dd\N) sea . (68) 

Thus, in the scalar channel there is a contribution of disconnected terms 



Pd£;,n=l-muu\0}(N\uu\N) sea (69) 



corresponding to the approximation, expressed by Eq. (|6*7j) . Of course, there is a similar 
expression for P^ n 

P!C,n = I ((0\uu\0)(N\dd\N) sea + (0\dd\0)(N\uu\N) sea ) . (70) 

4 Interference terms 

We now turn to the situation when one of qT x q operators acts on the constituent quark 
while the other one acts on the pion. In the one-loop approximation of the PCQM this 
contribution corresponds to the Feynman diagrams illustrated by Fig.3. 
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4.1 Contact interference 



The four-quark condensate can connect the pions of the nucleon self-energy loop with the 
quarks composing the nucleon. The contribution can be presented as 

WIT* T y \N) = V f (^ol g j|^^)^lg rX gl^ , )( 7r lg ry gl 7r )^" / ' 7r l ij 'j|^o) ( 71 \ 
[ W qq ql ' Q j^ n J (E -k 10 -E n + is)(E -k 2 o-E n , + ie) 1 } 

xA^k^A^h) 



(2tt)H (2tt)H 

where ki,k 2 are the four-momenta of the pions, A n (k) = l/(k 2 — M 2 + is) is the pion 
propagator. Recall that we include intermediate quark states with n = n' = only. The 
state vectors \<p n ) compose the complete set of the quark states with the energy E n , index 
corresponds to the ground state and Hi denotes the quark-pion interaction (f4*2|) . The 
summation is carried out over the quarks Q which compose the nucleon. 

In the quark language this means that the four-quark condensate can connect the 
pions with the intermediate quark of the self-energy loop or with another quark. These 
contributions are shown in Figs.3a,b. The corresponding exchange diagrams are shown 
in Figs. 3c, d. These expectation values contain the matrix elements (7c\qT x q\Tr) and 
(Q\qT Y q\Q). The former has a nonvanishing value in the scalar case only. The latter 
matrix element survives in the scalar and vector cases only, for a unpolarized nucleon. 
Thus, only the expectation values J ss and J sv obtain nonzero values. 

The connections of the pion n a with the intermediate state quark J xy (shown in 
Fig. 3a) and with another quark K XY (shown in Fig. 3b) are tied by the relation 

KY = -21^ (72) 

for a fixed quark flavor. This relation can be obtained by comparing the results of the 
integration over the pion energy in the loops of the diagrams shown in Figs.3a,c. 

One can write 

Ia Y = -(K a \qT x q\n a ) J d 3 zF a (z)ij(z)T Y ^(z)F ii (z). (73) 

with 

F a (z) = — J d 3 x^(x)ij B T a ^(x)S(x)D 7V (x - z) (74) 
where ^ is the SU(2) doublet of the light quarks. In Eq. (|74^l 

DM = j (T5) 

47T X 

is the three-dimensional pion propagator with \x = M n , while 

S(x) =M + cx 2 (76) 
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is the scalar field with the parameters JT] 

M = ^m' C= 2R3- (77) 

The parameters /3 and R already occurred in the definition of the valence quark wave 
function of Eq.(23). Only terms with the scalar structure T x provide a nonzero value. 
Also the integral in the rhs of Eq. (|73*j) does not turn to zero for the scalar and vector 
structures T Y only. 

The total contribution of such interference terms to the expectation values of the 
operators qq ■ qY Y q {Y = S, V) can be expressed by the contribution of the tt U loop 
to the self-energy diagram of U quark - see Eq.(|75j) 

J SY = -§ ■ (1 + f " - (78) 

with the factor (1 + y) taking into account the exchange diagram shown in Fig. 3c - 
Eq . (J44j) . while the coefficient | is the weight of the color-asymmetric state. The factor 
n takes into account the charge dependence of the ttQQ vertices and the number of the 
corresponding diagrams. We find n = 20 for the (mm) 2 condensate in the proton, turning 
to n = 7 for the neutron. It is n = 27 for the uudd condensate. For the scalar-vector 
condensate ddu^^u we have n = 20 for the proton and n = 7 for the neutron. 

Details of the calculation for the value Iq Y are given in Appendix C. Here we present 
the result. By expressing the pion matrix element by the vacuum one - see Eq. (JlOJ) . we 
obtain 

I SY = CjA Y (0\qq\0) (79) 

with 

A Y = fp\t)e-^ Y (t), (80) 

while 

ni\ - ,, , 

2 J w 2 - (3 2 

Here we use the standard notation erf(x) = J* e~ y2 dy and 

v s (t) = l-f3H\ V v (t) = l + f3H\ (82) 

being caused by the matrix elements ipT s ip = ipip and ^r v ip = ipjoip. For the coefficient 
Cj we get 

The interference terms provide for the scalar condensates 

J^' uu = 6.2 • 10~ 2 4 J^ uu = 2.2 ■ 10~ 2 4 J^ ud = J^ ud = 8.4 • 10" 2 ^. (84) 



m = 4erf(t)-te-t 2 -^-e-* 
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For the operator ddwy^u we finally obtain 

J^' du = 9.6 • 1(T 2 ^ and J* v > du = 3.4 • 1(T 2 ^ (85) 

where the lower indices p, n represent the proton and neutron, J^ua^ = Jm\o^o in the 
nucleon rest frame. 

Note that the insertion of the four-quark operator can lead to a charge-exchange 
pion-quark interaction between the points of emission and absorption of the pion by the 
constituent quark. This mechanism is also described by the diagram of Figs.3a-d and 
provides a contribution to the expectation value uT x ddT Y u. The charge-exchange matrix 
element (7r°|<i7 M 'u|7r + ) which is related to the vector part of the weak decay amplitude 
7r + — > 7i°e + i> e , has a non-zero value. However, the contribution is suppressed by an addi- 
tional small factor rrig/M^ when compared to the expectation values (ir a \qq\7r a ) . When 
calculated in the approach described in Sect. 3 the matrix elements (ir°\d''y p u\'n' + ) and 
(7r°|<iu|7r + ) vanish. Nonvanishing values are provided by corrections of the relative order 
M m to the PCAC relation expressed by Eq. (j4T?j) . 



4.2 Vertex interference 

Another type of interference term, illustrated by Fig.3e,f, is due to the PCAC relation 

(0\qi P l 5 r a q\^ a (k)) = iV2F T k p (86) 

where the pseudovector current connects the pion and vacuum states. The equations 
of motion lead to similar relations for the matrix elements of the pseudoscalar operator 
between the pion and vacuum states. In particular 

{0\u^d\n ) = ■ (87) 

m u + m d 

where q 2 denotes the square of the four-momentum of the pion. If one of the matrices 
acting on Lorentz indices, i.e. T Y , has a pseudovector or a pseudoscalar structure, there 
is a nonvanishing matrix element 

(<Po\qT Z T a q ■ gr y r a g|0o,vr a ) = (0 o |gr z r a g|0 o )(O|gr y r Q g|7r a ) (88) 

for any matrix T z and 

(<f) \qT Z T a q\(f) } = J d 3 x^(x)T Z T a ^(x). (89) 

The sum over color is carried out in both matrix elements in the rhs of Eq. (J88|) . 

Replacing the PCQM vertex in the nucleon self-energy by the vertex defined by Eq. 
we obtain (see Figs.3e,f) 
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Q n V E — k — E n + ie 

(O|gr^r a g |7r a )(0o|gr y T a g|0 n )(0 ra ;7r a |g f |0 o ) \ rf 4 fc 

The rhs of Eq. does not turn to zero only when both r x and T y are either pseu- 
dovector or pseudoscalar matrices. These cases must be treated separately. We shall use 
the standard PCQM approach [TQJ 120] where the sum over n in Eq.(|9"U]) is restricted to 
the quark ground state. 



4.2.1 Pseudovector case 

We start with the case where both matrices T x and T Y in Eq. (|9()|) are the pseudovector 
ones. The manifestation of the vertex interference in the self-energy diagrams is described 
by a certain tensor l pa with vanishing time components - see Appendix B, Eq.(Bll). Thus 
the integration over the energy k in the rhs of Eq. (j90j) can be carried out in the same 
way as for the self-energy PCQM diagrams [20] • We present the contribution as 

V = £J£> (91) 

a 

'F a (z)^(z) lal5 T a ^(z)(0\q(z) lpl5 r a q(z)\n a ) 

+F%z)^{z) 1 ^r a ^{z)^ a \q{z) lpl ,r a q{zm_ . 

The two terms correspond to the manifestation of the mechanism in the two vertices of 
the one-loop self-energy diagram - see Fig.3e. The functions F a (z) are determined by 
Eq.(J711). 

Since the operator qr°q does not change flavor, the neutral pions it contribute to the 
pseudovector condensates (uT A u) 2 and uT A udT A d only. The charged pions ^ provide 
contributions to the expectation values of the operators u^ p ^dd^ a ^u. Thus the charged 
pions give contributions to all the structures uY x udY Y d with the weights being determined 
by the Fierz transform. 

The tensor structure 

I2* = C?(g pa -Wf) (92) 

is determined by setting the time components to zero. The coefficients Cf = ^I pa g pa can 
be expressed through the contributions Eg (Fig.3e) of each quark to the total self-energy 
of the nucleon. 



r 



d A z 
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We start with the interference in the tt + UD vertex. Putting T x = 70-75, T Y = 7 P 7s in 
Eq.(|S5)l and projecting it on the quark states treated in momentum space, we obtain 

r d 3 k' - - - 
(U\d lplb u\-K + )(U\u lal sd\D)g(>° = V2F W J $(k')j p k p ij 5 iP(k' - k) (93) 

= 2y/2F n j ——^(k')i l5 S(k)^(k' - k) 

with S(k) standing for the scalar effective field, while k! denotes the momentum of the 
quark. The last equality is due to the PCQM equation of motion. The rhs of Eq. (|9*3^1 is 
2F% times the PCQM quark-pion vertex - Eq.p2|) 

Being substituted to Eq.flSH) for J ±1 , Eq.® provides the value Cf = f-F^Sj with 
Eg denoting the contributions of the single quarks in the self energies T, ± . One can also 
obtain that Cj = |FJEq. To show this, note that 

uTu = ^r(r° + I)q (94) 

(with / standing for the 2x2 unit matrix) for any 4x4 matrix V acting on Lorentz 
indices. Since (0|gr/g|7r°) = 0, one finds 

±=(U\u lplh u\A = (95) 

with the further procedure as in the charged case. 
All the contributions can be expressed by the value 

C? = -i^E , (96) 

which can be presented as 

°' 6n 2 (1 + 1/32)2 Ri ]„ "V y'+^f 1 ' 

In the chiral limit /i 2 = we obtain 

r°- ^ 1 = 1(32 + ^ 1 - 47 - 10 ' 3 - 1 14 1(TV (98) 

while for // = M n we have 

C\ = -1.0 • 10" 2 ^. (99) 

This provides 

2 10 

Jtr = 3 • + y) ■ (100) 
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where the factor (1 + also includes the exchange diagram shown in Fig.3f, - Eq. 
while the factor | is the weight of the color asymmetric state. Thus we finally have 

j^r = -1.4 ■ lo-wfc, - (ioi) 

Turning to condensates of quarks with different flavors, we set for the contribution of 
the neutral pions to the pseudovector structure 



C = (~n u - n d )(l + H) • | . tf^jj, = 4.2 • KT 2 ^ - ^ 



The charged pions provide a direct contribution to the expectation values of the oper- 
ators ujp^ddjcrj^u = J2 c ,g plhd c d 9 ^ a lhU g with c and g standing for the color indices. 
Their contribution to the expectation values of the operators u a T x u a 'd b T Y d b \5 aa ,5 w - 
8ab>fiba'), which we are looking for, is determined by the Fierz transform 

<'4 = -k E V % dT A u 5 a , b - i £ T%\l, b dT A \ K u. (103) 
In our case the pseudovector term with the diagonal color structure contributes only to 
uT x udT Y d = --•-• J 7p75 u uT x Yl^ Y d + ... (104) 

Here the summation over colors is carried out, providing the factor — ~. The dots denote 
the terms which do not contribute. 

Following the previous analysis we must separate the pseudovector component 70-75 in 
the operator r X 7 p 7 5 r y in the rhs of Eq. ([104|) . The interference terms can be expressed 
by the tensor 

1% = 2(1 + H) . ( nu + n d )I Q pa . (105) 

The tensor is obtained by the summation of the rhs of Eq. (|91j) over the charged 
pion states and over the constituent quarks of the nucleon and by inclusion of the exchange 
terms. The coefficient Cj is given by Eq. (j96j) . The contributions are 



J SM = -\g pa i% J Ps ' ud = -J SM , (106) 



jV,ud = K g jC _ 1 jC 
° fj,v o U pp u i^^ 1 pa piSiivilpo 1 pa i 



jA,ud _ jA,ud 1 r r jC _ \ jC 
' J pv °0pu ' ^UppUi/a 1 pa QiJpviJpa 1 pa 

T T,ud 1 / \ T C 

J uv,a3 = -^VKVuvlpVaplajLa- 
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4.2.2 Pseudoscalar case 

Now we consider the pseudoscalar case, i.e. T x = T Y = 75 in Eq. (j90j) . If the charged 
pions are exchanged the matrix elements of the quark operators between the vacuum and 
the pion states are given by Eq. (j87|) which respects the isospin symmetry. However, the 
contribution of the neutral pion exchange contains the matrix elements 

which depend on the quark masses m u ^ separately. This breaks for m u 7^ rrid the isospin 
symmetry explicitly. 

After the integration over ko in the rhs of Eq. (j90)l (see Appendix D) we can present 
the contribution in a form similar to the pseudovector case -Eq. (|5T|) 

i = J2i a ( 108 ) 

a 

Sz[F a (z)^(z) l5 T a i;(zm^)l^q(z)\K a ) 

+F a (z)${z)j 5 i;(z)(7r a \q(z) l5 T a q( Z )\0)} 

where we must put q 2 = m 2 in the matrix elements determined by Eqs. (|87J) . Q108|) . The 
expectation values can be expressed by the term I + of Eq. (jl08j) corresponding to the 7r + 
meson 

2M 2 

' d 6 xd 6 yS(x)ip(x)^ 5 ilj(x)D n (x - y)if>(y)jBip(y) (109) 



m u + m d 

with D w and S defined by Eqs. (J75)) . (|75jl . Numerically we get 

J+ = 7.0 • 10~ 2 4 (110) 

Proceeding in the same way as in the pseudovector case we find for the contribu- 
tions of the interference terms containing the neutral ir° mesons to (M75M) 2 and u^ud^d 
condensates 

O If) 

J^ s ' uu = -{1 + ^)nj u (111) 

and 

9 10 

jPsM = + _ ){n J d + n J u) (n2) 

with 

I + ■ (m u + m d ) 



Using m u = 4 MeV, = 7 MeV we find 

I u = 4.8 • 10" 2 4, I d = 2.8 • 10~ 2 4 



;il3) 
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The charged n mesons contribute to the expectation values of the operators u^dd^u 
providing thus the contributions to all basic structures defined by Eq.(j3J). They contain 
the factor 

Ic = (1 + + n d )I + = O.Ue 3 (114) 

being 

jS,ud = Ij C) jPs,ud = jS,ud + jPsM^ ( U5 ) 

6 

T V,ud _ _1 f jA,ud _ £ f 

jT,ud _ 2 ~ 

In more sophisticated models of the pions |2Sj the quarks obtain large effective masses. 
Thus the terms I a will become much smaller. 

4.2.3 Mixed case 

If one of the matrices in the rhs of Eq.()88|) is a pseudoscalar (75) while the other one 
describes the pseudovector (7 p 7s) we find contributions to the condensates with mixed 
Dirac structures qY x qqY Y q. If T x = 7 P 75, T Y = 75 the expectation value turns to zero 
when we neglect the possible intermediate state excitations of the constituent quarks - 
see Eq.(Bll). The terms with T Y = / y p / -f5,r x = 75 provide nonzero values. When we 
focus on the expectation value ddu / -f fl u (which is ddu^ou in the rest frame of the nucleon) 
among the mixed condensates, we must calculate the expectation value of the operators 
U'j 5 ddj 'j 5 u and u^o^dd^u. Contrary to the pseudoscalar case T x = T Y = 75, such 
terms do not contain the large factor M n /(m u + ma) ^ 12 - see Eq. (llU9|) . providing thus 
a minor contribution ~ 10~ 3 4 

4.3 Total contribution of the interference 

Now we can present the total contribution of the interference terms. For the quarks of 
the same flavor they are presented in Eqs.(jBU), P01|> - (|1H|) and 

= 0.064 J n uu = 0.024 (116) 

J Ps ' uu = 0.07n u e 3 , J^ u = -0.014n u (^ - ^)e 3 (117) 

turning to zero for the other structures. Recall that n u stands for a number of u valence 
quarks in a nucleon. For the quarks of different flavors Eqs. (j84J) . (|106|) . ()115|) provide for 
the proton 

jSM = 22£ 3 j jPsM = _ .284 (H8) 
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Jll ud = (-0.05^ + OM^el, J^ d = O.Ug^el 

JtSrt = (0-25 0.08^ 

Of course, the values presented by Eq. (jll8|) coincide for the proton and neutron except 
for the pseudoscalar case where 

Jp SM ~ Jn s,ud = 0.034 (119) 

Also the value for the (u^u) 2 condensate for the proton differs from the value (d^d) 2 for 
the neutron by 

(Jo S ' UU )p ~ (Jo S,d X = 0.064 (120) 

These characteristics obtain nonzero values due to the explicit dependence on the cur- 
rent quark masses. As we noted earlier, these effects would be much smaller if more 
sophisticated models for the pions are used 

Finally, for the scalar- vector condensates we obtain by using Eq. (}85|) 
jsv,du = 9 6 . 1(T 2 ^4 J^' du = 3.4 ■ 10^4 



5 The values of the four-quark condensates 

Now we sum the partial contributions obtained in the previous Sections. We present 
the results in units of the characteristic scale = — (0\qq\0) (see Eq.(|33j)\ The partial 
contributions to the expectation values defined by Eq.Q are due to the constituent quarks 
(denoted by C and shown in Fig.l), to the pion cloud (denoted by P and shown by Fig.2) 
and to the interference terms (denoted by J and shown in Fig. 3). This is expressed by 
Eq.pjl. We do not present the values of the parameters which are negligibly small in our 
scale. 

A. Scalar channel. Recall that in the scalar case there are specific "disconnected terms" 
in which one of the products qq acts on the QCD vacuum. Such terms emerge from 
the contributions of the pion cloud. Thus in the scalar case the values of the four- 
quark condensates can be presented as the sum of the "disconnected terms" Udis and the 
"internal terms" Ui nt 

U p (n) = Udis]p(n) + Ui n t-p(n) (121) 

with the indices p, n denoting proton or neutron. For the other structures the " discon- 
nected terms" vanish. 

We start by presenting the results for the disconnected terms. Using Eqs. (j63|) and 
(foTj) we obtain 

U S d r = P S d C U = -3-83^ (122) 

for both proton and neutron. 
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Consider now the "internal" contributions. For the scalar case they are expressed by 

Eqs.flHl), (EB, flUD, (ESI), (ED, (HUH) 



ttS,uu /~iS,uu I tjS,uu i jS,uu ri i i 3 /-i qo\ 

^intjp ~~ u p r int;p * J p ~ U.llfcg, l izo J 

ttS,uu /~iS,uu _|_ pS,uu | jS,uu pi nq-3 

u int;n ~ n ">" r int;n * J n ~ U.ZOfc Q , 

with the partial contribution 

qf<- = 0.084 c? w = 0.014 P s,uu = P s,uu = _ 25£ 3 ; (124) 

jV« = 0.064 J^ uu = 0.024 
The total values of U^ u (IT2TJ) are the sums of Ea. (1X221 and 

jjs,uu = _ 3Me 3^ u^uu = _ 4 _o5 e 3. (125) 

In this case the sea-quarks mainly contribute, the rest coming from the direct action of 
the four-quark operator on the constituent quarks and from the interference terms. 

Considering the mixed-flavor condensates uudd, we obtain for the disconnected terms 
presented by Eq. ffTUj) 

Ug* = PS? = -3.06eg (126) 

for both proton and neutron. The internal terms are expressed by Eqs. (j28)l . (jSHJ), ({EH}, 
and Ea. fTTHjl : 

c sm = QA6e ^ p SM = _ QM£ ^ jSM = 0>22e 3 ? ( 127 ) 

U^f = C SM + P%t d + J SM = -0.134 



^ = ^ = uSM + v sm = _ 319£l (12g) 
There are no "disconnected terms" in the other channels, thus U = Ut nt 



B. Scalar- vector channel. For this case there is no contribution if all the four quarks 
belong to the sea. They contribute through the interference determined by Eq. (j85|) while 
the constituent quark contribution is given by Eq.(jBH), (f3T)J) 

C sv ' du = 0.18e 3 (129) 

for both proton and neutron. The interference terms are ()85j) 

jSV,du = q ^q^ jSV,du = o 03£ 3 

Thus Eqs.(j3IJ), (jHOj) and (|83|) provide for the mixed scalar- vector condensate ddu^ou 

jjSV,du = C SV,du + jSV,du = 0.28£q, (130) 
v SV,du = C SV,du + jSV,du = q 2l£ 3 
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C. Pseudoscalar channel. For the pseudoscalar case we find by using Eqs. (j29|) . (jSJ, 

dsi, (Ei, dm 

^Ps.uu = ^Ps,™ + pPs,ud + jPs,uu = _ X g l£ 3 ; ( 131 ) 
^Ps.uu = ^Ps,™ + pPs,ud + jP S ,«n = _! g 6£ 3 

The partial values are 

Cf-™ = -0.024 pr*'™ = pP»>™ = -2.034, 

J^ s ' uu = Q.Uel, J£ s > uu = 0.074 

The numerical values are determined mostly by the contribution p Ps > uu of the sea quarks. 
In the case of the condensate uj^udj^d we obtain from Eqs. (j2l?j) . (J3~3~j) . (JIT]) . (ESI), (E3) 
and (dUI) 

C Ps ' ud = -0.03^, pPsM = 0-51^, (132) 
jPsM = _ 2 8el J? s > ud = -0.31^ 
composing, following Eq.(jlJ) 

jjPsM = c Ps,ud + pPs,W + jPs,ud = 0>20e 3 ? ( 133 ) 
jjPs.ud £jPs,ud _|_ pPs,ud _|_ jPs,ud g ]^£-3 

The difference of the values Up S,ud and U^ s,ud is caused by the explicit dependence on 
the current quark masses - see Eq. ()119|) . 



For the vector, axial and tensor structures we present the values of the coefficients 
a v ( A < T ) anc l b v ( A > T ). Recall that we introduced a notation where the partial contributions 
of the valence and sea quarks and of the interference terms are denoted by the lower 
indices C, P and J (see the text below Eq.fjSJ)). The second lower index labels the specific 
nucleon. Thus, ap p denotes the contribution of the sea quarks to the parameter a v of 
the proton, etc. The notations a5 n ) and &^ n ), labeling the vector (V), axial (A) and 
tensor (T) kept for the total contributions to these parameters for the proton 

(neutron). We omit the lower index if the values coincide for both nucleons. Note that in 
all the channels the sea-quarks do not contribute to the parameter 6 pn , i.e. 

&P,p(n) = (134) 

for all composition of flavors. Recall also that the interference does not contribute to the 
expectation values of the operator of the same flavors in the vector and tensor cases - see 
Sect.4. 

D. Vector channel. By using Eqs.(|30"|). (EEUi (O) . (IMJ) we obtain 

a v d ™ = -0.01^, = 4;r = -0.514, (135) 
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qV,uu _ a< j^ _)_ cipp — 0.524 dn' uu — dp n — 0.5l£o 

while 

h V,uu = b V,uu = Q 13£ 3 ) = b V,uu = q Q2£ 3 (1 g 6) 

For the mixed-flavor condensate we find IjBTty . (Pjl . pHjl . (jfiSj) . (fTT£l) 



r/ 



^ = _ .02eg, 4^ = 0.514 a? d = -0.054 (137) 



a^ d = a% ud + + = 0.444 
which are the same for the proton and neutron, as well as the parameters 

bl' ud = 0.254 h Y d = 0-044 (138) 
b VM = h VM + b VM = Q 20£ 3 

E. Pseudovector channel. Here we find by using Eqs. flSTJ), (gT]), fH, (EH), (fTTTIl 

a c 'p U = —0.034 ap'p U = a P '™ u = 0.5l4 (139) 

a% uu = -0.034 <r = -0.01eg, 
a A >™ = apT + a A P C + <AT = 0.45e 3 



p ~ u C,p u P,p u J,p ~ u.<±ue , 

-4,"" _ „A™ I A,uu , A,uu _ n r n 3 
n _ a C,n "T «P,n "I" a J,n ~ U.OUfcg 



and 



= 0.034 6jT = 0.034 = 0.014 (140) 

b A,uu = b A,uu + 6 ^«u = g £ 3 ; 
= b A,uu = Q Ql£ 3 

For the expectation value of the operator uY A udY A d we get with Eqs. (j31|) . ()33|). ()41|) . 

(El, (El, (PSD 

a^ = -0.06^, ^ = -0.514 = 0.144 (141) 

a A M = a A,ud + a A,„d + a A„d = _ 43£ 3 

and 

b A,ud = b A,ud = q Q6£ 3 ( M2 ) 

F. Tensor channel. Using Eqs.flH2J), (@TJ), (ESI), (El) we obtain 

T.uu n r\ a 3 T.uu T.uu no 3 /-1 /io\ 

«c, P = 0.044 a P ; p = a P ; n = -1.02eg, (143) 

„T,uu _ T,uu , T,uu _ n QOc-3 
a p — CLq p + CLp p — — U.yo^Q, 
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a n' UU — a C,n + a P,n ~ 1-02£q 

while 

b T p uu = 6g u = -0.054 (144) 
For the mixed-flavor operator flUHJ), (jUJ), flU, llTTB|) 

u £ ud = 0.074 ap ud = -0.51ejj, a T y ud = 0.254 (145) 

a T,«d = TM + a T,ud + a T,ud = _ Q19e 3 

while 

b T d ud = -0.104 b T / ud = -0.084 (146) 
b TM = b T,ud + b Tud = _ 0A8 3 



The final results of this Section are presented in a compact form in Tables 1,2 (keeping 
the values larger 0.1 in modulus). The numbers are given in the units of eg = 1.47- 10~ 2 
GeV 3 , see Eq.(J35J). The values U x > hh are 

U^ hh = (N\ : qf^T x q^' ■ q^ b T x q^' : \N){5 m ,8 w - 5 ab ,5 a , b ) 

with N = p,n - see Eqs.JTJ, ©, ®, (UK- 



Table 1 



X 


TjX,UU 
P 


TjX,UU 

n 


TjX,ud 


TjX,ud 
n 


s 


-3.9 


-4.1 


-3.2 


-3.2 


sv 






0.3 


0.3 


Ps 


-1.9 


-2.0 


0.2 


0.2 



Table 2 



X 


n X,uu 
Up 


I.X, uu 
U P 


n X,uu 
a n 


UX,uu 
n 


n X,ud 

a v 


hX,ud 
U P 


n X ,ud 
a n 


UX,ud 
n 


V 


-0.5 


0.1 


-0.5 





0.4 


0.3 


0.4 


0.3 


A 


0.5 


0.1 


0.5 





-0.4 


0.1 


-0.4 


0.1 


T 


-1.0 


-0.1 


-1.0 





-0.2 


-0.2 


-0.2 


-0.2 



6 Summary 

We calculated the expectation values of the four-quark QCD operators qY x qqY Y q in 
nucleons for all basic Lorentz structures and for compositions of the light quark flavors. 
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We employed previously derived results of the perturbative chiral quark model (PCQM) 
which treats the nucleon as a system of three valence quarks surrounded by a pion cloud. 
We approximate the averaging of the product of operators over the valence quark by the 
matrix elements of the constituent quark operators over the PCQM constituent quarks. 
We present the expectation values of the operators acting on the sea quarks by the ex- 
pectation values of QCD operators in pions. The intensity of the pion field is determined 
by the PCQM model result. 

The expectation values of the scalar and pseudoscalar operators are Lorentz scalars. In 
the other channels they have a more complicated tensor structure being determined by the 
two parameters a v( - A ' T ^ and b v ^ A ' T ^ -Eqs.dHJ), ©• For the quark operators with the same 
flavor, e.g. qT x qqT x q, the scalar and pseudoscalar condensates, as well as the parameters 
a v ( A > T ) for the other structures are dominated by the contribution of the sea quarks. The 
averaging of four [/-quark operators over the valence quarks in the neutron provide zero 
values in the lowest order of PCQM. This occurs because the operators uu should act on 
different quarks while there is only one £/-quark in the neutron. In the case of the proton 
both constituent quark and interference terms provide minor corrections of the order of 
several percent to the main contribution of the sea quarks. In the contrary, the sea quarks 
do not contribute to the coefficients b v ^ A ' T \ In the vector and tensor channels the values 
b v ' T for the proton are determined by the contribution of the constituent quarks. 

In the case of the mixed-flavor condensate uT x udT x d the role of the vertex interference 
increases due to the large combinatorial factor. These terms become as important as the 
sea-quark terms in most of the channels. The parameters b v ^ A ' T ^ are determined mostly 
by the contributions of the constituent quarks. 

The contributions of the sea quarks are expressed by the expectation values of the four- 
quark operators in pions. Latter values are in turn expressed by the expectation values 
of the four-quark operators in vacuum jT^j. The specific numerical values are obtained 
by using the vacuum factorization approximations [T3]. Thus the contribution of the sea 
quarks is expressed by the well known vacuum expectation value (0|gg|0). 

In the case of the scalar-vector condensate ddu^QU there is no contribution coming 
from the pions only. Averaging over the neutron is dominated by the contribution of the 
constituent quarks. The interference and the constituent quark terms are of the same 
order of magnitude in the proton. 

We can draw some conclusions on the chiral properties of the expectation values which 
we study in the present paper. The contribution of the sea-quarks has the same explicit 
dependence on the pion mass M n as the contribution of the sea-quarks to the expectation 
value (N\qq\N). The latter expectation value, which is proportional to M% times the pion- 
nucleon cr-term, is known to depend strongly on M n . On the contrary, our interference 
terms exhibit only a weak dependence on M n . The valence quark contribution does not 
contain an explicit dependence on M n . 

Due to the explicit dependence of the vertex interference terms on the quark current 
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masses, we have the isotopic symmetry breaking effects in the pseudoscalar channel. The 
absolute magnitude of this effect is numerically small with several units of the value 10 _2 £o 
for our scale Eq. The effect is much smaller if the quarks, composing pions are assumed 
to have the constituent (but not current) masses [23] . 

In the special case of the scalar condensate the expectation values are dominated by 
"disconnected terms" in which one of the quark operators acts "inside" the nucleon while 
the other one acts on the QCD vacuum. This contribution comes from the sea quarks, 
reflecting the pion structure [T5] . 

Note that a nucleon expectation value is the excess of the density of the quark operator 
products over the vacuum density, integrated over the volume of the nucleon 

(N\qT x qqT Y q\N) = (N\ J d 3 x[q(x)T x q(x)q(x)T Y q(x) - (0\qT x qqT Y q\0)}\N) . (147) 

Of course, the first term in the rhs of Eq. (|147j) is positive. However, the whole rhs of 
Eq. (|147|) can be negative. This is why some of the expectation values run negative. 

In an earlier calculation jH] the scalar expectation value {N\(uu + dd) 2 \N) was deter- 
mined in the framework of the Nambu-Jona-Lasinio model [7j under certain additional 
assumptions. Actually, the expectation values of the color-singlet four-quark operators 
q a q a q b q b have been obtained in [B] . Thus, to compare to the results of [6^ we must extend 
our analysis to such operators as well. 

In the paper jB] the expectation value is presented as the composition of the contribu- 
tion of the constituent quarks Ad and of a and it mesons, A a and A n . Our contribution 
of the constituent quarks appears to be several times smaller than the value of Ajj. The 
large discrepancy is not surprising, since the conception of the constituent quarks is quite 
different in the two models. The meson contribution A a + A n of 6J could be compared 
with the "internal" sea-quark contribution of the present model, containing the expecta- 
tion value, which is presented by the second term of the rhs of Eq. (jllj) . The corresponding 
contribution P (with the "hat" sign labeling the color singlet operator) can be obtained 

by using the formula obtained in [TSj. The result P = §J]§^|& should be compared to the 

sum Aa+AnOf®. We find P = 1.53s 3 , = 2.3- 1(T 2 GeV 3 while A a + A T = 3.6- 1CT 2 GeV 3 . 
The total values in the two models are A a + A v + Ap in [B| and the sum U = P + C + J 
in our approach. We obtain C = \C\ J — |j where the additional term J dominates 
in the sum C + J. The NJL value is A a + A w + A D = 6.4 ■ 1(T 2 GeV 3 while we obtain 
U = 2.5£q = 3.7 ■ 1CT 2 GeV 3 . The results provided by the two approaches differ by a 
factor of about 1.7. One of the possible reasons for the discrepancy is that some of the 
contributions have not been accounted for in both approaches. 



Two of us (E.G.D. and V.A.S.) are grateful for the hospitality of University of Tubingen 
during their visit. The work was supported by the DFG grants: 438/RUS 113/595/0-1, 
FA67/25-1, GRK683 and by the RFBR grant 03-02-17724. 
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7 Appendix A 



Here we show how the contributions to the four-quark expectation values, obtained in 
the paper, manifest themselves with the help of the PCQM formalism. As an example 
we consider the operator uuuu averaged over the proton. In the framework of the PCQM 
the nucleon is a system of three constituent quarks, where the bare three-quark state is 
renormalized by nN interactions. Thus, the physical proton state \N) is expressed as 

\N) = Texp(-i : f° dtH^t) :)|0 O ) 

where |0o) is the state of three valence quarks and Hj is the renormalized Hamiltonian 
of the interaction between the constituent quark and the pions which includes the coun- 
terterms. 

The expectation value {N\(uu) 2 \N) can then be written as 

(N\(uu) 2 \N) = Z 2 (<p \(uu) 2 \<p Q ) 
+2(0 o |M 2 |0o)(0o|^;i0o,7r)(0o,7r|^;i0o) 

+ (0 O |#;|0O, 7T> (00, 7T|(^) 2 |0O, 7T)(0o, A^M*) ■ (Al) 

Here Z = 1 + |^ is the renormalization factor, while S is the sum of the self-energies of 
the constituent quarks with energy E. 

In the next step we present each pair of the operators uu as the sum of operators 
acting on the valence and the sea quarks 

uu = (uu) v + (uu) s . (A2) 



Thus 

((j) \(uu) s \(f)o) = 0, (7r|(uu) v |7r) = 
and Eq.(Al) takes the form 

(N\(uu) 2 \N) = (1 + 2— )(0o|M 2 |0o> (A3) 

+2(0o|(m) 2 |0o)(0o|^/|0o,7r)(0o,7r|^/|0o) 
+ ((fro\H I \(j)o,7c)((fr \(uu) 2 v \(f) )((f) ,7r\H I \(j) ) 

+ (0o I Hi 1 0o , ?r) (vr | (uu) 2 S \ n) (0 O , ^ I Hj \ <f> Q ) 

+ 2 ( o | Hi 1 00 , 7T ) ( 0o | ( UU ) v 1 0o ) ( 7T I ( UU ) 8 \ TV ) ( , VT | \ ) . 



Since we include the quark-pion interactions to lowest order, the renormalization effects 

as 



are taken into account in the first term of the rhs of Eq.(A3) only. We put Z 2 = 1 + 2 dl 
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The rhs of Eq.(A3) can be simplified due to some cancellations. The second term 
in the rhs describes the self-energy insertions. These contributions are canceled by the 
counterterms of the PCQM Lagrangian JU]. Another cancellation occurs between the 
third term and the part of the first term 

2— (0 o |W'|0o) + (0o|(^/|0o,vr)(0o|W'|0o)(0o,vr|^|0o,7r) =0. (AA) 

This can be obtained in a straightforward way. The last equation is a rather standard 
cancellation of the radiative correction by the renormalization factor. Note that the two 
terms in the rhs of Eq.(A4) do not cancel totally for the operators (uu) 2 averaged over 
the neutron - see Sect. 2. 

Thus, Eq.(Al) takes the form 

(N\(uu) 2 \N) = <M^l0o> 

+(0o|^/|0o,7r)(7ri(ww)^|7r}(0o,7r| J ff/|0o) (A5) 
+2((p \H I \(p ,7r)((f) \(uu) v \(f)o)(n\(uu) s \n) (</> , 7r|i?/|0 o ) • 

Now we can identify the terms in the rhs of Eq.(A5). The first term corresponds to 
the contributions of the constituent quarks. The second term describes the contribution 
of the sea quarks shown in Fig.2. The third term presents the "interference" effects with 
one of the uu pairs coming from pions while another one comes from the constituent 
quark. The latter can be the same as the one in the matrix element of the interaction Hj 
or the other one. These terms are shown in Fig.3a-d. A cancellation similar to Eq.(A4) 
takes place for all the operators (uT x u) 2 in the proton, although in the general case the 
operator depends on the spin variables. However, the spin dependence manifests itself 
through the operator with / and II denoting the two quarks. Since the color 

wave function is asymmetric, the two quarks compose the spin-symmetric state being at 
the same space point. Thus, the two-quark spin wave function \x I)H ) is the eigenfunction 
of the operator (a 1 a 11 ) with (a 1 '^ II )\x 1 ' 11 ) — \x I,H )- Hence, the four-quark expectation 
value can be separated as a factor and the cancellation takes place as well as in the scalar 
case. Similar analysis can be carried out for the operators of the general form qT x qqT Y q. 

In the special case of the axial and pseudoscalar operators, there can be the interference 
effects in the first order of the ttQ interaction. This happens because the matrix elements 
(0|gr x g|7r) have nonzero values in these cases. Thus, the operators qT x qqT x q determine 
a nQQ vertex (Q\qT x qqT x q\Q , n) . This causes the "vertex interference" contributions 

(N\qT x qqT x q\N) mtrf = (0 o |#i|0o, (0 o , Ti\qT x qqT x qlfo) (AG) 

+ (<Po\qT x qqT x q\<j )o , 7r)(0 o , 7t|# 7 |0 o > 
with X labeling an axial or pseudoscalar. Such terms are shown in Figs.3e,f. 
In the rhs of Eq.(A6) we have 

(0o,vr|gr x ggr x g|0 o ) = (n\qT x q\0) (fo\qT x q\<p ) (A7) 
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(0o|gr x ggr x g|0 o ,vr) = (o|gr x g |vr)(0 o |gr x g|0o). 

We assume that the matrix elements of the QCD operators (qq) 2 over O are approx- 
imated by the matrix elements of the renormalized PCQM constituent quark operators, 

(<h\ml\<M = (4>o\(Q r Ql 2 \4>o) (AS) 

in the first terms of the rhs of Eqs.(A3) and (A5). The renormalization ^Hj means that 
the shape of the constituent quark wave function is modified by the influence of the pion 
cloud. Also, we approximate the matrix element (<f>o\(qq)v\<f>o) — (4 > o\QQ\4 > o) i n the third 
term of the rhs of Eq.(A5). 



8 Appendix B 

Except for the scalar case, the matrix element between the two-quark states depends on 
the spin orientation, containing the factor with / and II denoting the two quarks. 

Since the color wave function is asymmetric, the two quarks compose the spin-symmetric 
state when being at the same space-point. Thus we must put 

( x I >»\(a I a II )\ X I > II ) = l (Bl) 

for the value of the spin operator (a 1 a 11 ) averaged over the spin two-quark wave function 
X I,n of the quarks / and II. 

For the scalar case we find immediately 

with T(x) defined by Eq. (jl9j) . while g(x) = e~ 2x2 ^ R2 N 4 . This provides 

C& Utt = ^ 2 (l-|/3 2 + ^) (53) 
for the proton, with M defined by Eq . (|27p . For the pseudoscalar case we get 

T[x) = -4f3 2 g(xy- ^ 1 (S4) 

leading to 

For the vector and pseudovector structures we can find in the rest frame of the nucleon 



a V c iA) = -\C$ A % (B6) 
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with i and j being the space indices, corresponding to the four-dimensional indices ji and 
v. A direct calculation provides for the vector case 

= ^ 2 9(x)^ 2 ■ l^a 11 ) - Jo* 1 ), (57) 



leading to 

To determine the coefficient b^, we calculate the time components 



a v c =N\-\n (58) 



and, since Coo = + 6^, we find 

b v c =^(1 + ^ + ^(3% (BIO) 
o lb 

For the pseudovector case notice that the time components turn to zero. We introduce 
the notation 

K = ( v&x ) 

for the bispinor entering the wave function - Eq. (j23j) . We obtain 

^7o75^ = 0. (511) 

Thus 

a£ + b£ = 0. (512) 

As to the value of a A , it can be calculated by using Eq.(B6). In the pseudovector case we 
get 

R lilh K = Vi + (3 2 K J R ; (513) 
By using the properties of the Pauli matrices one finds 

(<jx)ai(ox) = 2xi(ax) — x 2 Oi. (514) 

Thus 



= 9 (,) (o / + + p W^-?^ {m5) 

leading to 



^=A/- 2 (-i + i/3 2 -lA (510) 
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Finally, in the tensor case we find for the space components 

2 (ax)a k (ax) 

KCJijK = E ijk (a k - (3 — ) (517) 

and the function / can be obtained by using Eq.(B14). For the space-time components 
we have 

Ra 0j K = -2(3^ (518) 
R 



and 



F<>3,ok = ^PV(x)^6 jk (519) 



with the further procedure described in the main text. 



9 Appendix C 

In order to calculate the value Iq Y introduced by Eq.(|7S|). we must calculate the function 
F°(z) - Eq . (I74|) . For the w-quark it takes the form 

F°(z) = / d 3 x X *(ax) X S(x)D n (x - z)$ 2 (x), (CI) 



F n R 

changing the sign for the <i-quark. We present the pion propagator (J7HJ) as 

/J3t. p ik{x—z) p —ixa 

The factor e~ l( - xa \a = 0) is introduced in order to simplify the calculations by expressing 

xD n (x — z) = iV a D w (x — z). (C3) 
We obtain by doing the integral over x 

7T 3 / 2 N 2 BR 4 

F(z) = -- f^x^tVMAT^z) + BT 2 (z)), (b = 0) (C4) 

with A = M + §aR 2 , B = -jcR 4 , while c and R are determined by Eas. l|2H]) . (j77|) and 

r d 3 k p -iHz~b)-\k 2 R 2 
r rP>h f,-ik(z-b)-^k 2 R 2 
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We can evaluate the rhs of Eqs.(C5), (C6) by presenting 



1 ' " dae- a{k2+ ^\ (C7) 



k 2 + /i 2 Jo 



leading to 



CO 



T\{z) = J" da J _0_ e ^-S)-|^^--a^ + ^) 5 {C8) 

T 2 (z) = ^ 2 ±e^ - ^ Tl (z). (C9) 

The further calculations can be simplified by assuming the chiral limit fi 2 = 0. The 
integral in the rhs of Eq.(C8) is dominated by the values of z 2 close to |i? 2 . Thus, the 
integral over k 2 is determined by k 2 ~ while the integral over ct is dominated by 
a ~ f-R 2 - Hence, the factor afi 2 in the power of the exponent in the rhs of Eq.(C8) is 
about 0.12. Since Iq Y provides a small correction only, this makes the calculation of this 
value in the chiral limit fi 2 = reasonable. 

Calculation of the integrals over k and over a (by the substitution t = (\R 2 + a) 5 ) 
leads to Eqs.d?HJ) - JED of the text. 



10 Appendix D 

The integrals over the time component ko in the rhs of Eq. (j90|) take the form 

x= fdko q 2 1 (m) 

J 2ni q 2 -M 2 +ie E -k -E n + ie 1 ' 

with q 2 = k 2 , — k 2 . We can present X = X\ + X2 with 

Xi = Ml I ^ • = \ , (D2) 

W J 2iri q 2 -M 2 + ie E - k - E n + ie v ; 

X 2 =[^ -i . (D3) 

J 2ni E Q — k — E n + ie y 1 

The integral X2 can be expressed through the contribution of the pole in the upper half- 
plane of the complex variable ko- This corresponds to the negative-energy solutions of 
the Dirac equation. Such terms are neglected in framework of the PCQM. Hence, we put 
X = Xi, leading to Ea. (fTU3J) . 
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11 Figure captions 



Fig.l. The contribution of the valence quarks to the expectation values of the four-quark 
operators. Solid lines denote the valence quarks, the dark squares denote the four-quark 
operator. 

Fig.2. The contribution of the sea quarks to the expectation values of the four-quark 
operators. The dashed lines represent the pions. Other notations are the same as in 
Fig.l. The dark circles denote the vertices of the pion-quark interaction. 

Fig.3. The contribution of the interference term to the expectation values of the four-quark 
operators. The "contact interference" is illustrated by Figs.3a-d. The "vertex interfer- 
ence" is shown in Figs.3e,f. The permutated diagrams are not shown. The notations are 
the same as in Figs. 1,2. 
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